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Density functional theory methods are applied to investigate the properties of the new superconductor (3- 
YbAlIU and its polymorph a-YbAlB4. We utilize the generalized gradient approximation + Hubbard U 
(GGA+U) approach with spin-orbit(SO) coupling to approximate the effects of the strong correlations due to 
the open 4/ shell of Yb. We examine closely the differences in crystal bonding and symmetry of /3-YbAlB4 
and a-YbAlB4. The in-plane bonding structure amongst the dominant itinerant electrons in the boron sheets 
is shown to differ significantly. Our calculations indicate that, in both polymorphs, the localized 4f electrons 
hybridize strongly with the conduction sea when compared to the related materials YbRh 2 Si2 and YbB 2 . Com- 
paring /3-YbAlB4 to the electronic structure of related crystal structures indicates a key role of the 7-member 
boron coordination of the Yb ion in /3-YbAlB4 in producing its enhanced Kondo scale and superconductivity. 
The Kondo scale is shown to depend strongly on the angle between the B neighbors and the Yb ion, relative to 
the x — y plane, which relates some of the physical behavior to structural characteristics. 

PACS numbers: 71.27.+a, 74.70.Tx, 71.15.Mb 



I. Introduction 

Heavy electron systems have provided key insights into emer- 
gent quantum mechanical behavior in correlated materials. 
Such systems may be tuned by the application of pressure 
or a magnetic field through a quantum critical poinP ^ In 
this region they may exhibit non-Fermi liquid propertie d 1 * 2 ^ 3 ! 4 ! 
and often superconductivity 56 . Quantum criticality ha s been 
widely reported in Ce, U and Yb systemP ™ 10 1 11 1 12 1 . How- 
ever, while associated heavy-fermio n superco nductivity has 
been reported in Ce and U compound J 7 * 8 ! 9 ! 1 ° l 1 -^ 2 4 examples in 
Yb systems have been lacking. The discovery of superconduc- 
tivity alongside non-Fermi liquid behavior 13 in the clean Yb 
based heavy fermion (HF) compound /3-YbAlB 4 has there- 
fore excited much interest. Thus there exists a clear impetus 
to determine the electronic structure of /3-YbAlB 4 and its re- 
lationship to other materials. 

The superconductivity and pronounced non-Fermi liquid 
behavior of /3-YbAlB 4 occur at ambient pressure 13 and zero 
applied magnetic field. At zero field the resistivity p ex T 3 / 2 
(T is the absolute temperature) and there is a logarithmic di- 
vergence in the specific heat. The introduction of an external 
magnetic field tunes the system away from quantum criticality 
and towards conventional Fermi liquid behavior p ex T 2 . Re- 
cent quantum oscillation studies demonstrate heavy fermion 
behavior in the material 14 . Furthermore, the Curie-like sus- 
ceptibility at high temperatures suggests the presence of lo- 
calized moments that are pivotal to the formation of the heavy 
Fermi liquid. 

/3-YbAlB 4 is similar to its polymorph, a-YbAlB 4 , in ex- 
hibiting Curie-Weiss type behavior in the susceptibility at 
high temperature and in the formation of a heavy fermion 
ground state at low temperatures. /3-YbAlB 4 differs from its 
polymorph, a-YbAlB 4 , firstly in that superconductivity has 
only been found in /3-YbAlB 4 (below 80mK). Furthermore, 
measurements of the low temperature specific heat 15 show a 
larger linear specific heat coefficient saturates at above 150 
mol -1 K _1 with the application of small magnetic fieldP 



in /3-YbAlB 4 compared to approximately 130 mol 1 K 1 a- 
YbAlBjSl. This may point to stronger correlation effects in /3- 
YbAlB 4 that contribute to its superconductivity. Electron spin 
resonance measurements 16 also suggest that /3-YbAlB 4 pos- 
sesses a unique signature of both local moment and conduc- 
tion electron behavior that has been attributed to its quantum 
criticality and is not observed in a-YbAlB 4 . Also, while fj- 
YbAlB 4 shows features of quantum criticality and non-Fermi 
liquid behavior at low temperatures and ambient pressure, a- 
YbAlB 4 behaves as a conventional Fermi liquid 17 . 

In this investigation we consider the electronic and mag- 
netic properties that may drive the intriguing behavior of 
both /3-YbAlB 4 and its polymorph a-YbAlB 4 . We utilize 
density functional methods with DFT+U to model the system. 
Importantly, these calculations differ from the itinerant / 
electron calculation favored in a recent quantum oscillation 
study of the low temperature Fermi surfaced However, a 
DFT+U approach allows us to consider the driving interac- 
tions behind the behavior of the system, in particular in a 
Kondo physics framework. We also consider the relationship 
between these materials and related 4f 13 materials such as 
YbRli2Si2 and YbB2. In particular we consider the coupling 
of the 4f moments to the conduction sea and the proximity of 
/3-YbAlB 4 to magnetic instabilities. 

II. Crystal Structure 

In this investigation we use the crystal structure extracted 
experimentally by Macaluso et all 15 for both /3-YbAlB 4 and 
a-YbAlB 4 . /3-YbAlB 4 crystallizes in the ThMoB 4 structure 
with space group Cmmm and a-YbAlB 4 crystallizes in the 
YCrB 4 structure with space group Pbam. Both space groups 
are orthorhombic. The lattice parameters are shown in Ta- 
ble [I] We also show the crystal structures of both polymorphs 
in Fig. [T] We have doubled the conventional cells of each 
structure along the a-axis to aid visualization of the bonding 
networks. The lower symmetry structure of a-YbAlB 4 in the 
lower panel may be most easily compared to /3-YbAlB 4 via 
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to its nearest boron atom. In /3-YbAlB 4 these bonds lie along 
a single axis while there are two different orientations in a- 
YbAlB 4 . Although this difference is unlikely to influence the 
orbital state of an individual Yb Atom, it may be influential 
upon the magnetic coupling between the Yb 4/ magnetic mo- 
ments. The alteration in symmetry is then likely to influence 
the magnetic state of the system and in this way be important 
to understanding the properties of these materials. 

In this study we also consider the properties of two related 
crystal structures that both have 6-member coordination of the 
Yb ion. The first is a hypothetical structure, which we re- 
fer to as 7-YbAlB4, and has the same c-axis lattice param- 
eter and cell volume as /3-YbAlB 4 , but a honeycomb boron 
sublattice that leads to a regular 6-member coordination of 
the Yb and Al ions. The second structure we compare to is 
YbB2(YbYbB4) which exists in the same hexagonal symme- 
try as MgB 2 (P6/mmm) ls . The cell parameters for both of 
these structures are also shown in Table [I] Experimentally in 
YbB 2 the Yb f 13 orbitals exhibit a Curie- Weiss like suscep- 
tibility at high temperatures and the compound is thought to 
order antiferromagnetically 18 at T/v = 5.6 ± 0.2K. 



III. Electronic Structure Methods 



FIG. 1: (Color online) Crystal structure of (top) Cmmm f3- 
YbAlB4 and (bottom) Pbam a-YbAlB4 viewed along the c-axis The 
green rectangular outlines show the conventional unit cells and we 
have doubled the unit cell along the a-axis. The black rectangle in 
the lower figure shows an outline of a unit cell in a-YbAlB4 that 
may be most easily compared to the conventional cell of the upper 
/3-YbAlB 4 . 



the construction of the approximate unit cell outlined in the 
lower panel of Fig. [T] In comparing this new rotated cell in 
a-YbAlB 4 with the structure of /3-YbAlB 4 , the bonding net- 
works become similar. However, subtle differences in sym- 
metry mean that the transport, heat capacity, and the bonding 
in the two systems are different. 

Both structures are orthorhombic with relatively short c- 
axis lattice parameters which will bear importance for our 
electronic structure results. Both structures also contain a 
layer of interconnected 5 -member and 7 -member boron rings 
intercalated with Al and Yb atoms. Each Yb atom is coor- 
dinated above and below by 7-member boron rings in both 
structures. These seven member rings dominate the crystal 
field environment of the Yb atom and therefore impact its or- 
bital state. 

Crystallographically, the two polymorphs are primarily dis- 
tinguished by the manner in which the 5-member and 7- 
member rings of the boron plane are interconnected. Both the 
5-member and 7-member rings are irregular polygons. This ir- 
regularity is far more pronounced in the a structure. Therefore 
the intercalant Yb and Al atoms do not reside equally close to 
all boron atoms in the nearest boron ring in the two structures. 
For example, in Fig. [T] the bonds from the Al show the link 



The electronic structures determined in this investigation are 
based on the all-electron approach to density functional the- 
ory using full-potentials in WIEN2lP. The Brillouin zone 
integration was achieved by the tetrahedron method typically 
using al8 x 18 x 30 k-point mesh. The radii of the muffin 
tins were set at 1.6ao for B, 2.11ao for Al and 2.15ao for Yb. 
We utilised RK max = 8.0. 

We have calculated the electronic structure using the GGA- 
PBE 20 correlation functional and with +U on the Yb sites for 
the 4f orbitals. We employ U = 8 eV and J = 1 eV and spin- 
orbit coupling in a second variational approach. The 'around 
mean field' double counting correction in WIEN2k was used. 

In this section we compare and contrast the electronic struc- 
ture of /3-YbAlB 4 and a-YbAlB 4 with 7-YbAlB 4 and YbB 2 . 
We begin by showing the calculated density of states for /?- 
YbAlB 4 and a-YbAlB 4 in Fig. [2] In all cases the large peaks 
in the DOS are dominated by the Yb 4f contributions. The 
overall structure of the DOS for /3-YbAlB 4 and a-YbAlB 4 are 
similar. For both, a large manifold of states due to the filled/ 
orbitals lies between -7 eV and -3 eV below the Fermi level. 
The states between -2 eV and +2 eV about the Fermi level 
are dominated by itinerant states that have little charge within 
the muffin tins. Visualization of the charge density shows that 
most of the itinerant charge is associated with the boron lay- 
ers. The empty 4/ 13 hole state sits just beyond 2 eV above 
the Fermi level. It is possible to obtain GGA+U solutions 
with a variety of different hole states (however, not all can be 
achieved). We find that the J z = —5/2 hole orbital always 
leads to the lower energy. 

The J z = —5/2 hole is consistent with the result pre- 
dicted by an approximate treatment of the crystal field sym- 
metry of the Yb site 21 . A low lying doublet of mj = ±5/2 
was predicted and shown to give a good fit to the measured 
magnetic susceptibility at high temperatures. Furthermore, 
our total energy calculations suggest that the energy scale for 
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TABLE I: Space group, cell parameters, atomic positions and site symmetries for a-, ft- and 7- YbAlB 4 as well as YbB 2 . 







a 


P 


7 


YbB 2 


Space group 




Pnam (No.55) 


Cmmm (No. 65) 


Pmmm (No. 47) 


P6/mmm (No. 191) 


Lattice 


a 


5.9220 


7.3080 


3.1346 


3.2561 


constants 


b 


11.4370 


9.3150 


5.4293 




(A) 


c 


3.5060 


3.4980 


3.4980 


3.7351 


Volume (A 3 /f.u.) 




59.55 


59.53 


59.53 


34.29 


Atomic 


Yb 


Ag (..m) 


4/ (m2m) 


la (mmm) 


la (6/mmm) 


positions 




(0.1294, 0.10543,0) 


(0, 0.30059, 0) 


(0, 0, 0) 


(0, 0, 0) 


(Wyckoff 


Al 


Ag (..m) 


Ag (2mm) 


If (mmm) 




position, 




(0.1387, 0.4096, 0) 


(0.1816, 0, 0) 


(1/2, 1/2, 0) 




site 


Bl 


Ah (..m) 


Ah (2mm) 


2p (m2m) 


2d (-6m2) 


symmetry, 




(0.2893,0.3126, 1/2) 


(0.1240, 1/2, 1/2) 


(1/2, 1/6, 1/2) 


(1/3, 2/3, 1/2) 


coordinates 


B2 


Ah (..m) 


8q (..m) 


2n (m2m) 




(x,y,z)) 




(0.2893,0.3126, 1/2) 


(0.2232, 0.1609, 1/2) 


(0, 1/3, 1/2) 






B3 


Ah (..m) 
(0.3840, 0.0468, 1/2) 


Aj (m2m) 
(0, 0.0920, 1/2) 








B4 


Ah (..m) 
(0.4740, 0.1943, 1/2) 
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FIG. 2: (Color online) Density of states for /3-YbAlB 4 (top) and 
a-YbAlB 4 (bottom). 



non-Ising flops from this state is at least 0.3eV (~ 348 IK) 
which implies that thermal flops in this system are likely to be 
Ising-like. Calculating the expectation value of the Ising-like 
moment gives 3.4/i# which is comparable to the measured 
values of the effective moment of 2.9 fi B 15 and 3.1/i#l In 
our calculations the order is ferromagnetic and the magneti- 
zation is directed along the c-axis which is the easy axis in 
the measured magnetization 15 . We note that the calculated 
specific heat coefficients of 1.21 and 1.54 mJ/mol K 2 for (3- 
YbAlB 4 and a-YbAlB 4 correspond to large mass enhance- 
ments ~100. This implies the presence of a low temperature 
heavy fermion ground state that is not represented by the bare 
DFT+U calculation. 

In Fig. [3] we show the DOS for 7-YbAlB 4 and YbB 2 . Both 
exhibit a similar localised Af 13 hole structure, but the hole 
state lies significantly further from the Fermi level than for 
both polymorphs of YbAlB 4 at approximately +3 eV. We 
find the entire manifold of /-states, both above and below the 
Fermi level, is shifted upwards and this is a result of the dif- 
fering crystal field produced by the coordination of 6-member 
boron rings as well as the effective volume about the Yb site. 
Interestingly YbB 2 seems to exhibit a very nearly gapped den- 
sity of states at « +1.75 eV. Together these indicate that the 
Aj hole may be less strongly linked by hybridization and spin- 
orbit effects with the itinerant states near the Fermi level. 

IV. Competing States and Quantum Criticality 

We digress from the electronic structure to discuss briefly 
some analysis relating to heavy fermion physics. In the Do- 
niach picture of heavy fermion systems two scales may be 
considered to dominate the physics: the Kondo scale Tk of 
local moment screening and the scale for the Af moments to 
magnetically order such as Trkky- By multiplying by the 
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FIG. 3: (Color online) Density of states for 7-YDAIB4 (top) and 
YbB 2 (bottom). 



density of states at the Fermi level these may be expressed in 
dimensionless form as: 



p(Ef)Tk oc exp 



p(E F )T RKKY oc ( 



where Jp(Ep) 



En 



En 



En 



(1) 



(2) 



(3) 



Here A and En are respectively the width and position with re- 
spect to the Fermi level E F of the 4/ level. Here J is the local 
moment - itinerant electron coupling as in the Kondo model. 
We note that this Doniach approximation does not account for 
the degeneracies of the 4f state in Yb 22 . The type of mag- 
netic orde r and di mensionality may have further influences on 
the system 2 MU. However, within the model we may gain 
insight into the drivers of the Kondo physics in these systems. 

Theoretical^ and experimental work using both doping 
and pressure tuning 17 indicate that /3-YbAlB4 is proximate 



FIG. 4: (Color online) Density of states for YbB 2 at the expanded 
unit cell volume: 46.68A 3 . 



to an antiferromagnetic instability. Experimentally, stoichio- 
metric YbB 2 is thought to order antiferromagneticallyS^ at 
Tat = 5.6 ± 0.2 K. Given the similar crystal structures of 
YbB 2 and /3-YbAlB4, we may compare these two compounds 
to determine which characteristics may promote the antiferro- 
magnetism in YbB 2 and the quantum criticality of /3-YbAlB4. 

In order to tune YbB 2 away from antiferromagnetism to- 
wards a quantum critical point we need to effect an enhance- 
ment of Tk over Trkky • Due to the exponential dependence 
of Tk this may be achieved by an increase in the quantity 
Jp(Ep) = We may, for example, tune this quantity by 
the applica tion of pressure, p. The application of pressure in 
Yb systemJ 2 ^ 7 -^ will increase A due to increased hybridiza- 
tion and also increases En since the valence fluctuation is be- 
tween Yb 2+ and Yb 3+ . Essentially, this is because the fluctu- 
ation involves the addition of an / electron and therefore the 
effective energy barrier to the process, £"0, is smaller when the 
volume about the Yb ion is larger. Experimentally, the behav- 
ior of En dominates 29 30 and we are tuned deeper into the anti- 
ferromagnetic state with increasing p. Therefore, we typically 
require an effective negative pressure, or volume increase, in 
order to tune from a magnetic state towards a quantum critical 
point. Experimentally, an effective negative pressure has so 
far not been achieved by chemical pressure in these systems, 
but we may simulate the effect by increasing the unit cell vol- 
ume in our electronic structure calculations. 

In Fig. [4] we show the density of states for the structure of 
YbB 2 at an expanded unit cell volume of 46. 68 A 3 . The vol- 
ume about the Yb ion is increased. As a result En falls and A 
also decreases compared to YbB 2 at its experimental volume 
in Fig. [3] In this structure the value of En is approximately 
the same as that in /3-YbAlB 4 . However, the width, A, in the 
expanded YbB 2 structure is much smaller than that found in 
/3-YbAlB4. As a result with such a volume tuning, A is likely 
to be small by the time we tune to a critical point. 

In Table jn] we compare the value of for several impor- 
tant structures. An important feature from the table is the abil- 
ity of /3-YbAlB4 and a-YbAlB4 to maintain a large hybridiza- 
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tion, A, even with a larger Yb-B distance compared to YbB 2 . 
This difference suggests that parameters beyond the volume 
about the Yb ion are important. The angle of inclination of 
the Yb ion to the boron ring may be pivotal to producing an 
enhanced (and Kondo scale) while maintaining a large A. 

To probe the importance of the angle, 0, subtended by the 
Yb ion to the boron ring we utilize our hypothetical 7-YbAlB 4 
structure. We do this by altering the aspect ratio of its unit cell 
to vary 0, while retaining a constant unit cell volume to create 
a set of structures that we denote as 7*-YbAlB4. In Fig. [5] we 
plot the value of A /Eq as a function of 0. We see that the peak 
value is centered about 52°, which is the angle subtended by 
the proposed ground state hole state \mj = ±5/2 >. Im- 
portantly, while the angle is altered we find little variation 
in Eq and it is the large changes in the hybridization, A, that 
drive the strong dependence of A/Eq on angle. Therefore, f3- 
YbAlB4 and a-YbAlB4 enhance their Kondo scale not only 
by increasing the effective Yb-B distance (to decrease Eq), but 
also approach the optimum angle for hybridization. It is the 
alternation of 5-member coordinated Al sites with 7-member 
coordinated Yb sites that allows /3-YbAlB 4 and a-YbAlB 4 to 
form this advantageous angle of inclination. In this sense 
we may consider /3-YbAlB4 and a-YbAlB4 to be symmetry 
tuned Kondo enhanced materials. In /3-YbAlB4 this Kondo 
enhancement drives the material to quantum criticality. 

Interesting future work with a fuller treatment of many- 
body effects with dynamical mean field theor)&^2] may at- 
tempt to demonstrate this relationship between crystal struc- 
ture and the heavy fermion ground state as has been shown 
in elemental Yb under the influence of pressure^. We also 
note that due to the nature of the DFT+U method, the absolute 
value of the parameter A / Eq may be sensitive to the value of 
U used in the calculations. However, the trend in A/ Eq as a 
function of the geometric change is consistent for different 
values of U. 

In this framework, we may also discuss why /3-YbAlB4 is 
the only discovered example of an Yb based heavy fermion 
superconductor, while many Ce based examples exist. We 
proceed by considering pressure tuning resulting in uniform 
volume contraction. The effect of pressure on will be dif- 
ferent for Ce and Yb compounds. For Ce compounds the ap- 
plication of pressure decreases Eq, but increases A. This is 
because the valence fluctuations are between Ce 3+ and Ce 4+ . 
This fluctuation involves the removal of an / electron and 
therefore the effective energy barrier, Eq, is small when the 
volume about the Yb ion is smaller. In contrast the application 
of pressure in Yb based compounds increases both Eq and A. 
In some sense these quantities are working against each other 
in terms of their Kondo physics in Yb compounds. For this 
reason at the quantum critical point, it is more likely that a 
narrow A will exist. Since A is indicative of the coupling to 
the itinerant electrons, then a lower tendency to superconduc- 
tivity may be found in Yb based systems. 

V. Band structure 

Critical to the operation of heavy fermion systems is the in- 
teraction of the localized /-electrons with the conduction sea. 
This interaction leads to the dramatic mass enhancements and 
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FIG. 5: (Color online) Plot of the quantity A/E versus the an- 
gle subtended by the Yb ion to the surrounding boron ring in 7*- 
YbAlB4. A/Eo is important to the determination of the Kondo scale 
and Table [n| indicates that the structure of /3-YbAlB4 lies nearest to 
the optimum. 



itinerant /-electrons in these systems via the Kondo effect. 
The RKKY interaction may also give rise to a significant scale 
for magnetic coupling that drives their magnetic properties. 

In Fig. [6] we show the band structure for /3-YbAlB4. We 
have elected to plot the dispersion along the symmetry di- 
rections for the reciprocal conventional cell, (r = (0 0), 
A = (1/2 0), B = (1/2 1/2 0), C = (0 1/2 0), 
Z = (0 1/2), Az = (1/2 1/2), Bz = (1/2 1/2 1/2) 
and Cz = (0 1/2 1/2), in Cartesian coordinates. For f3- 
YbAlB 4 there are two Yb atoms in the unit cell which lead 
to the appearance of two flat Yb bands at approximately +2.2 
eV above the Fermi level. The remaining 4/ states are lo- 
cated in a manifold between -7 eV and -3 eV below the Fermi 
level. With SO coupling these 4f bands are spin mixed and 
split into a 4/ 5 / 2 complex and a 4/ 7 / 2 complex separated by 
the spin-orbit splitting of roughly 1.5 eV. These manifolds are 
further split by the anisotropy of the Coulomb interactional. 
This along with a Hund's rule splitting due to the unoccupied 
4/-hole makes identifying the 4/ 5 / 2 and 4/ 7 / 2 states a chal- 
lenging task. 

Two important features of the band structure have been 
highlighted in Fig. [6] The red circle, centered at 1.5 eV at 
point C, indicates the location of a strongly hybridized /-hole 
band. In particular the archway shaped band shows strong 
/-character in a partial charges analysis. This /-character 
remains strong throughout its dispersion in the plane, and 
demonstrates a strong hybridization between the /-hole and 
the conduction states. In contrast such hybridization is al- 
most entirely absent in the related quantum critical compound 
YbRh 2 SijE(. 

In this band structure of /3-YbAlB4 the Yb moments are in 
a ferromagnetic configuration. This allows us to see the in- 
teraction of the Yb moment with the bands associated with 
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TABLE II: The value of the parameter . To estimate A we have used the full width at half maximum of the 4f hole peak in the density of 
states. 



Structure 


Yb-B Distance (A) 


Yb-B Angle 


A 


E 


A 

En 


/3-YbAlB 4 


2.70 


48° 


0.25 


2.2 


0.11 


a-YbAlB 4 


2.73 


48° 


0.23 


2.2 


0.10 


7-YbAlB 4 


2.516 


46° 


0.2 


2.8 


0.07 


YbB 2 


2.65 


45° 


0.2 


3.2 


0.06 




FIG. 6: (Color online) Band structure for /3-YbAlB4. The red cir- 
cle indicates the location of a strongly hybridized /-hole band. The 
green circle shows the location of strongly spin-split bands near the 
Fermi level (many splittings are not visible). Both indicate strong 
interaction of the /-hole with the itinerant electrons. 



itinerant electrons. The green circle, centered just above Ep 
at point B, shows the location of strongly spin-split bands near 
the Fermi level. Both indicate strong interaction of the local 
/-hole with the itinerant electrons. This has the benefit of pro- 
viding a measure of the degree of Kondo coupling of the Yb 
moment to the Fermi surfaces, because the exchange splitting 
of these bands reflects the coupling of the local moment to the 
itinerant bands. The exchange splitting of this band is approx- 
imately 100 meV, far larger than that found in YbRh 2 Si2 34 . 
Consistent with experiment, the T*, coherence temperature, 
in /3-YbAlBj^ is ~ 200K, almost an order of magnitude 
larger than that in 

YbRhsSiJ^M 

This large Kondo coupling 
may be critical in the behavior of these materials and is con- 
sistent with large estimates of T* from the measured specific 



heaP. As can be seen, however, the band crossing the Fermi 
level at Cz possesses very little exchange splitting. Therefore 
the Kondo coupling is highly anisotropic and likely to vary 
around the Fermi surface. This is consistent with the varying 
mass enhancements inferred from recent quantum oscillation 
studied. 

In Fig. [7] we show the band structure for a-YbAlB 4 . Again 
the Yb moments are in a ferromagnetic configuration. We 
note that the band structure has twice as many bands as j3- 
YbAlB 4 due to the fact that its primitive unit cell contains 4 
formula unit cells for a-YbAlB 4 , while for /3-YbAlB 4 there 
are only two. This leads to the appearance of increased com- 
plexity, yet the general structure of the two cases is similar. 
The position of the /-hole states at +2.2 eV above the Fermi 
level and the location of the 4/ 5 / 2 and 4/ 7 / 2 multiplets be- 
tween -6.5 eV and -3.5 eV below the Fermi level are similar. 
This similarity is due to the similar local symmetry of the Yb 
environment. The a-YbAlB 4 band structure also shows evi- 
dence for spin split bands near the Fermi level. An example 
is provided by the spin split band at T approximately +0.8 eV 
above the Fermi level. 

In Fig. [7] we show the band structure for YbB 2 . In this case 
there is only one formula unit per unit cell. Therefore there are 
fewer bands in the band structure and the Yb moment is nec- 
essarily in a ferromagnetic configuration. The band structure 
is highly dispersive in all directions. As expected from the 
density of states the /-hole level lies further from the Fermi 
level than in the cases of /3-YbAlB 4 and a-YbAlB 4 . Also, 
there is less evidence for the presence of spin split itinerant 
bands near the Fermi level. These features are indicative of a 
smaller Kondo scale in YbB 2 . 

VI. Charge Density Analysis 

In Fig. [9] we show contour plots of the valence charge den- 
sity in the boron plane for /3-YbAlB 4 , a-YbAlB 4 and YbB 2 . 
The dominant chemical bonding resides within this boron 
plane and is driven by the short boron to boron distances, all 
of which are less than 1.9 A. The boron bonding of these 
structures is very stable. For the high symmetry case of 
YbB 2 shown in Fig. [9jc) we see that the bonding between 
the boron sites in its 6-membered rings is the same for every 
B-B bond. This bonding is of a similar character to that in 
graphite. In contrast for /3-YbAlB 4 and a-YbAlB 4 all B-B 
bonds are not equivalent. The 5 -member and 7 -member rings 
determined by the crystallographic positions do not give rise 
to a dominant ring structure in the bonding charge density. 
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FIG. 7: (Color online) Band structure for a-YbAlB4. 




FIG. 8: (Color online) Band structure for YbB2 in the hexagonal 
Brillouin zone. 



Specific bonds within these 5 -member and 7 -member rings 
are stronger than others. 

As shown in Fig. |9ja) the dominant B-B bonds in f3- 
YbAlB4 form a 2D network of large loops each of which 
contain two formula units. The boron network is well con- 
nected in both the a-axis and b-axis directions. In Fig. |9jb) 
a-YbAlB 4 shows a different bonding network. The dominant 
B-B bonds form chains along the a-axis with double-boron 
branches. The difference in bonding anisotropy should also 
be reflected in anisotropy in transport behavior. As a result 
the bonding structure and distribution of the valence electrons 
can not be approximated as conforming to the ring structure 
as is the case in YbB 2 . 

The bond charge isocontours indicate that the bond popula- 
tion of the dominant bonds is approximately double that of the 
weak bonds. Interestingly, in both /3-YbAlB 4 and a-YbAlB 4 , 
the weak B-B bonds always lie between a Yb and an Al ion. 
There is a direct correlation between bond length and bond 
population. The weak bonds are all ^1.86 A in length, while 
all other B-B bonds are <1.77 A. Therefore the subtle sym- 
metry changes between these compounds are significant to the 
electronic structure. 

The different valence bond networks demonstrated to 
exist in /3-YbAlB 4 , a-YbAlB 4 and YbB 2 are likely then to 
affect the magnetic coupling between Yb moments. This 
may impact the RKKY interaction in these systems which 
may be the driving force behind magnetic instabilities. The 
uneven bond population surrounding the Yb ion also makes 
its environment anisotropic. This may be a key factor in 
producing the quantum criticality and superconductivity in 
/3-YbAlB 4 . 



VII. Conclusions 

We have obtained strong evidence for a large Kondo scale 
driven by strong interaction between the Yb /-state and the 
itinerant states in /3-YbAlB 4 and a-YbAlB 4 . This is evi- 
denced by strongly spin split bands and significant hybridiza- 
tion between the /-states and itinerant states near the Fermi 
level. 

We present a case for a key role of the angle subtended by 
the Yb ion to the boron ring in producing the large Kondo 
scale and possibly the superconductivity of /3-YbAlB 4 . The 
introduction of the 7 -member coordination of the Yb ion in 
/3-YbAlB 4 allows the angle of w 48° to be subtended to the 
boron ring. This angle brings the compound close to the opti- 
mal ~ 52°, for the maximization of A/Eq which promotes its 
large Kondo scale while still allowing for a large hybridization 
of the 4/-hole with the conduction states. The importance of 
the angle subtended by the Yb ion to the boron ring suggests 
a significant scope for the use of experiments involving pres- 
sure tuning via uniaxial pressure that may be able to replicate 
these conditions in other materials. The compound YbB 2 is a 
strong candidate for such measurements. The in-plane bond- 
ing structure amongst the dominant itinerant electrons in the 
boron sheets has been found to differ significantly between 
/3-YbAlB 4 and a-YbAlB 4 . 



(a) /3-YbAlB 4 



VIII. Acknowledgments 




(b) a-YbAlB 4 




(c) YbB 2 




FIG. 9: (Color online) Charge density contours in the boron plane 
on a linear scale for (a) /3-YbAlB 4 , (b) a-YbAlB 4 and (c) YbB 2 . 
For YbB 2 we have doubled the unit cell in the a and b directions to 
aid visualisation of the bonding. Note the different topology of the 
strong B-B bonds in the a and /3 types, both contrasting with the 
regular, hexagonal topology of YbB 2 . 
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